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Abstract 

We investigate the thermal conductivity in the armchair and zigzag M0S2 nanoribbons, by 
combining the non-equilibrium Green's function approach and the first-principles method. A strong 
chirality dependence is observed in the thermal conductivity. Particularly, the thermal conductivity 
for the armchair M0S2 nanoribbon is about 25% lower than that in the zigzag nanoribbon at room 
temperature. By calculating the Caroli transmission, we disclose the underlying mechanism for 
this strong chirality dependence to be the much less phonon transport channels in the armchair 
M0S2 nanoribbon in the frequency range of [150, 300] cm . 
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The two-dimensional single-layer M0S2 (SLMS) has drawn considerable interest in re- 
cent years, owning to its direct band gap and graphene-like layered honeycomb lattice 
structure.- - — The band gap in the SLMS is a great advantage over graphene, which has 
zero band gap without the help of strain or other gap-opening techniques. Thus, the SLMS 
becomes a promising alternative to graphene in many electronic fields. For this purpose, 
many efforts have been devoted to comparatively investigate various properties of the SLMS 
and the graphene. It was found that the interface between the SLMS and the gold metal is 
inefficient for electron injection; while the Ti metal can greatly improve the efficiency- The 
mechanical and electronic properties of the SLMS have been investigated by first-principles 
calculations, where the defects were found to play an important role.- Sanchez and Wirtz 
studied the phonon properties of the SLMS, using the first-principles calculations.- 

In a very recent experiment, Sahoo et.al measured the thermal conductivity in the SLMS 
through analyzing the laser power dependence of the Raman mode.- The obtained thermal 
conductivity is 52 Wm _1 K _1 at room temperature. The thermal conductivity is an impor- 
tant parameter for possible applications of the SLMS. High thermal conductivity is helpful 
to deliver heat away from the electronic device, which is good to prevent heating induced 
breakdown of the device. For thermoelectric applications, lower thermal conductivity can 
yield in higher figure of merit ZT.— ~— Hence, it is useful to provide an accurate prediction 
for the intrinsic thermal conductivity value, i.e the thermal conductivity in a pure SLMS 
without any phonon decaying mechanism. This intrinsic thermal conductivity can serve 
as an upper limit for the values from the experimental measurement or other theoretical 
calculations. 

In this letter, we comparatively study the thermal conductivity in both armchair and 
zigzag M0S2 nanoribbons (MSNRs), by using the ballistic non-equilibrium Green's function 
(NEGF) formalism and the first-principles calculations. We find that the thermal conduc- 
tivity is strongly dependent on the chirality. The room temperature thermal conductivity 
in the armchair MSNR is about 25% lower than that of the zigzag MSNR. The Caroli trans- 
mission function shows that the origin for the strong chirality dependence is because of the 
much less phonon transport channels in the armchair MSNR in the frequency range of [150, 
300] cm" 1 . 

In the ab initio calculation, we use the SIESTA package^ to optimize the structure of the 
SLMS. The local density approximation is applied to account for the exchange- correlation 
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function with Ceperley-Alder parametrization^ 3 . and the double-^ basis set orbital is adopted. 
During the conjugate- gradient optimization, the maximum force on each atom is smaller than 
0.01 eVA" 1 . A mesh cut off of 120 Ry is used. Periodic boundary condition is applied in the 
heat current direction and the in-plane transverse direction. The free boundary condition 
is applied to the out-of-plane direction by introducing sufficient vacuum space. There are 
two types of calculations in this work. In the calculation of the phonon specstrum, a small 
unit cell is used and an 8 x 1 x 1 Monkhorst-Pack k-grid is chosen for the sampling of the 
quasi-one-dimensional Brillouin zone. In the investigation of the thermal transport, a large 
unit cell is calculated and the Gamma point k sampling is adopted. 

The relaxed structures for armchair and zigzag MSNRs are shown in Fig. [TJ The Mo-S 
bond length is 2.414 A, which is between the value 2.382 A from Sanchez et.afc and 2.42 A 
from Ataca et.al— There are 36 (18) atoms in the smallest translation unit cell in the armchair 
(zigzag) MSNR. Fig. [2] compares the phonon dispersion in the armchair (left) and zigzag 
(middle) MSNRs. Two smallest translation unit cells are included for the armchair MSNR 
in the calculation of the phonon dispersion. For the zigzag MSNR, four smallest translation 
unit cells have been used. Hence there are 72 atoms involved in the calculation, resulting in 
216 phonon branches in both panels. The wave vector k is along the heat current direction 
and in the unit of 2ir/a, with a = 10.705 A or 12.364 A in the armchair or zigzag MSNR, 
respectively. An obvious band gap around 300 cm -1 is found in the phonon spectra for both 
armchair and zigzag MSNRs. Similar band gap has been found in some existing works.- 1 ^ 
The three acoustic phonon branches from the armchair and zigzag MSNRs are compared in 
the right panel. We note that the x axis here is in the unit of A" 1 , which is different from 
the left and middle panels in the same figure. It is obvious that the armchair MSNR has 
considerably higher phonon spectrum than the zigzag MSNR. It should be noted that the 
out-of-plane transverse branch should behave flexibly due to the two-dimensional nature of 
the SLMS. This flexure nature disappears in present calculation, due to the loss of the rigid 
rotational invariance symmetry in the ab initio calculation.— 

There are various options to investigate the thermal transport. Classical results can 
be obtained from the molecular dynamics simulation of the thermal transport, where the 
phonon-phonon scattering dominates the transport process. From the theoretical point 
of view, it is a big challenge to provide accurate prediction for the thermal conductivity, 
because the samples in the experiment always possess various unpreventable defects. Hence, 
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a more practical task is to provide an accurate (quantum) prediction for the upper limit 
of the thermal conductivity. For this purpose, we would rather apply the ballistic NEGF 
method.— 1 ^ It is based on quantum mechanics. The phonon-phonon scattering is ignored in 
the ballistic transport region, which is actually quite reasonable for low- dimensional nano- 
materials. For instance, in graphene, it has been found that the out-of-plane transverse 
phonon mode can transport almost ballistically even in a large piece of graphene.— The 
combination of the NEGF and the ab initio calculation can provide us an accurate upper 
limit for the thermal conductivity.— 

In the NEGF approach, the thermal conductance is calculated by the Landauer formula: 



where h is the Planck's constant. n(w,T) is the Bose-Einstein distribution function. The 
transmission T[u] is obtained from the Caroli formula: 



where G r is the retarded Green's function. G a = (G r ) is the advanced Green's function and 
Tl/r is the self-energy. These Green's functions can be calculated from the force constant 
matrix from the ab initio calculation.— 

Fig. [3] top panel shows the transmission functions for the armchair and zigzag MSNRs. 
These functions exhibit some regular steps, due to the absence of phonon-phonon scattering. 
In the ballistic transport, a is proportional to the cross-sectional area, since there are more 
channels available for heat delivery in thicker nanoribbons. However, the thermal conduc- 
tance from Eq. [T] does not depend on the length of the system, a can be used to get the 
thermal conductivity (re) of a MSNR with arbitrary length L: k = crL/s, where s is the 
cross-sectional area. We have assumed the thickness of the MSNR to be 6.033 A, which is 
the space between two adjacent layers in the bulk M0S2.- This thickness value is the same 
for both armchair and zigzag MSNRs, so its value does not affect our comparison for their 
thermal conductivity. It is quite obvious that the thermal conductivity does not depend on 
the cross section. It means that the thermal conductivity in the MSNR does not depend on 
its width, since the thickness is a constant. On the other hand, the thermal conductivity 
in ballistic region increases linearly with increasing length, which has been observed in the 
thermal conductivity of two-dimensional graphene.—^— In the experiment, the thermal 
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T[u>] = Tr (G r T L G a T R ), 



(2) 
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conductivity is measured for a MSNR sample of dimension around /im.- Hence, we will 
predict the thermal conductivity for a MSNR of length L = 1.0 /im. 

Fig- El bottom panel gives the thermal conductivity for both armchair and zigzag MSNRs 
of 1.0 jum in length. A strong chirality dependence is revealed in the thermal conductivity. 
The thermal conductivity in the armchair MSNR is smaller than that of the zigzag MSNR 
in the whole temperature range. The chirality-induced anisotropy is about 25% at room 
temperature. This anisotropy is much stronger than that in the graphene. In the ballistic 
region, the chirality-induced anisotropy is only about 1% in the single-layer graphene.— 

The top panel in the figure discloses the origin for this strong anisotropy. It is because the 
transmission function in the armchair MSNR is much smaller than that of the zigzag MSNR 
in the frequency range [150, 300] cm -1 . It means that the phonon transport channels are 
much less in this frequency domain for armchair MSNR, leading to much smaller thermal 
conductivity in the armchair MSNR, although the three acoustic velocities in the armchair 
MSNR are higher as shown in Fig.[2]right panel. Particularly, the room-temperature thermal 
conductivity is 673.6 and 841.1 Wm _1 K _1 in armchair and zigzag MSNRs. These values 
are about one third of the superior thermal conductivity value (around 3000 Wm _1 K _1 ) in 
the graphene,— because the overall phonon spectrum in the MSNR ([0, 500] cm" 1 ) is scaled 
by a factor of one third from the phonon spectrum in the graphene ([0, 1600] cm" 1 ). The 
experiment value is 52 Wm _1 K _1 , which is far bellow our quantum upper limit. It indicates 
strong phonon scattering or substrate-induced phonon leaking effects in the experimental 
samples, as suggested by the authors in Ref. |8. 

We should stress that out theoretical value serves as an upper limit for the thermal 
conductivity of the 1.0 /im MSNRs. If the experimental samples are of high quality, then 
the measured thermal conductivity should approach 673.6 or 841.1 Wm _1 K _1 from lower 
side. 

In conclusion, we have applied the ballistic NEGF approach to predict the upper ther- 
mal conductivity value in the armchair and zigzag MSNRs. The force constant matrix is 
calculated from the first-principles method. A strong chirality dependence is found in the 
thermal conductivity. More specially, the thermal conductivity in the armchair MSNR is 
about 25% smaller than the zigzag MSNR at room temperature. The Caroli transmission 
function discloses the origin for this strong anisotropy to be the much less transport channels 
in the armchair MSNR in the frequency range of [150, 300] cm -1 . 
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FIG. 1: (Color online) The relaxed structure of the armchair (top) and zigzag (bottom) MSNRs. 
Top and side views are shown for each system. Atoms are denoted by blue (Mo) and red (S). The 
translation unit cells are highlighted on the left end. 
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FIG. 2: (Color online) Phonon dispersion in armchair (left) and zigzag (middle) MSNRs. Insets 
show the unit cells used in the calculation of the phonon spectrum. The three acoustic phonon 
branches in armchair and zigzag systems are compared in the right panel. We note the band gap 
around 300 cm" 1 in the spectrum. 
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FIG. 3: (Color online) Top: the transmission functions for armchair and zigzag MSNRs. Note 
obvious less transmission channels in the armchair MSNR in [150, 300] cm^ 1 . Bottom: the thermal 
conductivity (k) obtained from thermal conductance (<r) through their relationship, k = aL/s, with 
cross-sectional area s = wh. The thickness h = 6.033 A, and the length L = 1 //m. a is calculated 
by the Landauer formula. 
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